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Recently, we addressed the problem of the correspondence between two different descriptions of the phenomenon of 
inelastic electronic transitions^between two distinct set of continuum electron states, through an intermediate resonant 
state coupled to an oscillatoiu. Even though the problem of interest there was the influence of vibration damping on 
the spectroscopy of adsorbates with the scanning tunneling microscope, it was shown that the formalism adopted, 
taken over from a description of resonant tunneling in semiconductor structuresB with the necessary modifications, 
was in fact reducible to a simpler resolvent method, in the no-damping limit, but with the restriction of just one 
vibration coordinate. Our intent at that point was to make a connection with a previous paperu on STM, but an 
ly-j . obvious possibility then was the eventual extension of that correspondence between the two approaches, from just 
' one vibration coordinate to a Brioullin zone of modes. The present paper addresses this issue, showing that indeed a 
simple resolvent model, without appealing to many-body field theory, seems to be able to describe inelastic tunneling 
OA ■ in small semiconductor structures, at least with the restrictions of a K temperature and an Einstein band of optical 
phonons. Even so, we believe it is of some interest to present this analysis, as it may provide a simpler means of 
assessing the essential physics of the situation and also as a starting point for further useful extensions. 

On the other hand, even though rthis model can be applied to many real distinct situations in physics and in 
chemistry, as pointed out by Gadzukul when drawing similar comparisons, we will be referring specifically to inelastic 
resonant electronic transport in semiconductors (for a recent exposition based on nonequilibrium Green functions, see 
Ref. ^|), so that our results can be compared directly with the more general, and more powerful, many-body approach 
used in Ref. |^, which employs a two-particle many-body Green's function and which was subsequently simplified to 
a singe-particle many-body Green's function treatments. Other methods have been published, such as those of Ref. 
^, for example. But with the present, simpler, method even the inclusion of more than one intermediate, resonant 
state should,not introduce major difficulties. The essentials of this method were used many years ago by Domcke and 
CederbaumEl, to describe inelastic electron scattering by molecules in the gas phase. 
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C<3 1 II. RESONANCE GREEN'S PHONON OPERATOR 

As usual, the total spinless Hamiltonian is a sum of three parts, the electron Hamiltonian, the phonon Hamiltonian 
and the electron-phonon interaction, H = H e i + H p h + Hint- 



H el = e°4c a + ^2 e l c i°l + ^2 ekC l Ck 

l k 



1 



k 



(Vakclck + V ka cic a ), 



(1) 



H ph = ^?i Wg (6tfe 9 + l/2), 



(2) 



and 



(3) 



In H e i the first term describes the localized resonant electron state, the second and third terms describe the electronic 
states in the left and right (isolated) leads, and the last two terms represent the hybridization between the localized 
electronic resonance and the two leads. H p h describes the harmonic phonon states and Hi nt is the linear electron- 
phonon interaction, whose strength is given by the coupling constant M q , acting only when an electron occupies the 
intermediate, resonant state. 

We proceed to compute the resonance Green's function from the corresponding resolvent operator G = (e — H) -1 , 
using the equations-of-motifln method, as in the original Anderson description of localized magnetic states associated 
with impurities in metalsau. In doing this, we assume the validity of the Born-Oppenheimer approximation and 
include the interaction Hamiltonian in the resonant electron state energy, writing it as a phonon operator 



e a = e a (b^) = el+J2M q (bl + b. q ), 



(4) 



thereby rendering the electron Hamiltonian H e i adiabatically dependent on the state of the lattice oscillators. 
From the equations of motion for the operators 



ih 



dc a 

at 



(5) 



we obtain, taking into account the commutation of electron and phonon operators under the adiabatic hypothesis, 
the set of equations 



ih 



dc a 
dt 

. dck 



(t) +J2 V akCk(t) +J2 V *'Mt) 



= efcC fc (i) + V ka c a {t) 

at 



it) 



dci 
dt 



eiCi{t) + Vi a c a (t). 



Fourier transforming to the energy domain, we get the corresponding algebraic system 

ec Q (e) = e a c a (e) + KfcCfe(e) + V a ici{e) 

k I 
ecfe(e) = e k c k (e) + V ka c a {e), Vfc, 
eci(e) = e/q(e) + y/ Q c a (e),VZ, 

which can be rewritten in matrix form 

el - | Vka e k 

Via 61 





(6) 



(7) 



(8) 



where is a column matrix, V^ k the row matrix transpose of Vak, eu and ei are diagonal matrices with elements 
e k t , £fe 2 j ' ' ' an d , £; 2 7 ' ' ■) respectively, 1 is the identity matrix of the appropriate dimension and arc null matrices 
of the appropriate dimension also. 

This matrix equality can be put in the form 



(el - F e i) c = 0, 



(9) 
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writing the previous matrices more compactly, in an obvious notation. Introducing now another diagonal phonon 
operator matrix 

F ph = \^hu q {b\b q + l/2)]l, (10) 
we define an elcctron-phonon resolvent operator 

G total (c) = [(c + ir?)1 _ (Fel + Fph)] -l 

= [(C + *77)l-F total ]-\ (11) 

C being the total energy parameter and r\ a positive infinitesimal. Then [(( + ir])l — F total ] ® G total (£) = 1 which, 
expanded, means 

(C + iri )-e a -H ph l -Vj k -V£ 

-V ka « + irj)l - e k - H ph l 

-V ta (C + - ei - H ph l 

G aa G ak G al \ / 1 \ 

Gka G kk G k , = 1 . (12) 
Gu Gi k Gn / V 1 / 

Introducing a common index m denoting both k and I, we have the system of equations 

[(C + if)) - Co - H ph }G aa - ^2 V am G ma = 1 

m—k.l 

-VmaGaa + [(C + "?) ~ Cm ~ H ph ]G ma = 0. (13) 

In the second equalitiy we solve for G ma (where m — k,l) and substitute it in the first, obtaining the Green's phonon 
operator for the localized electronic resonance 

G aa (C)= {[(( + iv)-e a -H ph ] 

- ]T \V am \ 2 [{Q + ir ] )-e m -H ph ]- 1 }- 1 (14) 

m—k,l 

as a function of the total energy parameter, £. At zero K temperature, the initial energy in the transition process, equal 
to the total energy, will comprise the incoming electron energy, e ki plus the ground state lattice energy, 1/2^ %us q . 
With this last quantity fixed, we use as the energy parameter the incoming electron energy and write 

G aa (e ki ) = {[{e ki + 1/2 £ hu q + if?) - e° - £ M q (b\ + 6_ 9 ) - £ hu q (b\b q + 1/2)] 

9 9 9 

-^|K™| 2 [e fei + 1/2^7^ + ir}- ^^ g (6t6 g + l/2)]- 1 }. (15) 

m q q 

However, since in general (x + irj)^ 1 = P(x^ 1 ) — in5(x), P designating the principal part, we get 

[e* 4 + iV ~ em ~ ^2 ^g^M" 1 

9 

= P[(e kt -€ m -^2hu q b\b q )- 1 ] 

9 

-inS(e ki -f m -^ hw q b\b q ) (16) 

9 

and defining an electronic shift A and width T phonon operators by 

A(e fci )=P \Vam\ 2 {e ki -e m -Y J H b i)~ 1 

m=k,l q 

r(e fei ) = 2^^ \V am \ 2 5{e ki - £ro - ^ fc&t^), (17) 

m—k,l q 
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we can rewrite the Green's operator for the resonant state as 

G aa {e ki ) = {[e ki -e° a -J2 M q (b\ + &_,) - £ huo q b\b q ] 
q q 

-A(e ki )+ l -F(e ki )}-\ (18) 

At this point, and for our present purposes, we neglect the operator nature of both A and T as well as their dependence 
on e, obtaining a simplified Green's operator 

G aa (e ki ) = [e* -e«-^M,(&t +&_,)- ^fo^&t&J- 1 (19) 
with a complex resonance energy e a = e° + A — iT/2, corresponding to a lorentzian profile. 



III. THE TRANSITION AMPLITUDE 



Having obtained the resonance Green's operator, we now proceed to compute the transition amplitude, from the 
T matrix operator, T = V + VGV, between an initial electron-phonon state \ki', n q = 0, Vg) with zero phonons and a 
final state \lf\ n,i, n 2 , ■ ■ ■ , njsf), with n q excited phonons in mode q, in a total of N lattice vibration modes. Since we 
assume no direct transitions between electron k states and / states, we have 



(lf\T\k t ) = (lf\(J2 Vla4c a )G aa (J2VakC f a C k )\k t 



(20) 



still a phonon operator, to be inserted between the vibrational ground state of the lattice and all possible excited 
states: 



(If, nx,n 2 , ■ ■ ■ , n N \T\ku 0i, 2 , • • • , On) = (ni,n 2 , ■ • ■ , n N \(l f \T\k i )\0 1 ,0 2 , • • • , On) 



Vi }i aVa,ki(ni,n2, ■ ■ ■ ,n N \ 



|o>, 



(21) 



where we wrote the lattice initial, ground state as |0). 

Next, we proceed to diagonalize this last denominator by appealing to the operator U (as in Refs. ^ and [l0]) 



U = exp 



(22) 



and inserting U 1 U — 1 in the matrix element (3.2) above, after noting its effect on the phonon operators, b Q and b q , 



(23) 



results that assume symmetry of the phonon bands, specifically that M q = M_ q and tu q = u>— q . We have also employed 
the well known operator theorem that exp{A}B exp{—A} = B + [A,B] if [A, B] is a c-number. Then, the matrix 
element in (3.2) can be written 



(ni,n 2 , • • • ,n N \U 1 U 



C/ _1 t7|0> 



(24) 
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and diagonalizes to 



(ni,n 2) - • -,n N \U 1 



J q u q u q 



U\0). 



(25) 



^- e ~°+E^-E^ 

<? y 9 

We proceed to calculate this matrix element. First we note that, since M q = M_ g and uj q = u_ q , we can write 

[/|O) = e -£*^- 6<!) |0i,()2,---,0jv> 

= e «=ST l6 «i 9l) |0i> • e »^2 {b «2 b ^)| ^.... 

M "N f( ,t _ b ) 

■■■e "=57^ 9jvJ |0jv> (26) 
from which, taking advantage of the operator rule exp{A + B} = cxp{A} exp{_B} exp{ — 1/2 [A, £?]}, we obtain 

'^(*i-*.)i 0ff > 

= e a U«J V A ^ | m ) ; 

z * / i yit\ 



g 1", 



m q — 



(27) 



for each vibration mode q. We then obtain the product of similar terms for all N phonon modes. Operating on that 
product with the diagonalized denominator phonon operator, we have 

-i -l 



i ~ e a ~ ^2 TtoJqb\b q 



U\0 1 ,0 2 ,---,0 N ) 



= e -5E,(4) v E ••• E m< 

mi— rri2— niN— 

(Mi/nu>i) mi (M 2 /nw 2 ) m2 • • • {M N /hu N ) mN 



\/m\\m-£. ■ ■ ■ mjv! 



\m 1 ,m 2 , ■ ■ -,m N ). 



(28) 



It remains to calculate the action of the operator J7 1 on the intermediate vibrational state ket |mi, m 2 , • • • , mpf): 



|mi,m 2 ,---,TOjv) 



(29) 



a product of terms of the type 



• r E E (-)' 



(M q /hb> q ) r <+ a < y/m q \{m q -r q + s q )\ 



=0 r„=0 



(m q -r q y. 



m q -r q + s q ) 



(30) 



Performing the internal product with the corresponding bra (n q \, we use the orthogonality condition, (n q \m q — r q +s q ) = 
S(n q — m q + r q — s q ), and get, for the q th mode, a contribution (including the pre-factor in Eq.(3.8) above) 



aU " J E - — 7— r Nl^K) 



m Q — 



2 /M D 



h-l/2 



E(-r«E(-) r ' 



m Q — 



(n q \){(M q /hu; q ) 2 } r * 
r q \(m q - r q )\(n q + r q - m q )V 



(31) 
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However, this last summation is a generalized Laguerre polynomial! 

TJ._I 



13L 



r„=0 



r g !(m g - r q )\(n q + r q - m q )\ 

= L%- m °[(M q /hu q f} 

and the whole transition amplitude matrix becomes 
(lf,ni,n 2 ,- ■ ■,7ijv|T|& i ;0;L ) 02 ) ' ■ • , Ojv) 



(32) 



V lf , a V a , ki e 



huii 



oo oo 



E E • • ■ E (-) mi ^r mi [(M 1 /fw.) !t ]{-) m °L%- m ° [{M 2 /hu 2 f\ ■ . ■ (-)«<■» L ™- m » 



{M N /huj N ) 



mi— mo—O rriN—0 



e, + E 



Ml 

TlLU n 



(mihuJi + m 2 hu 2 + ■ ■ • + toat^wat) 



(33) 



isary to distinguish the two cases, n q < m q and m q < n q , if we take into consideration the symmetry 



It is not n 

propertiescflc3 of the generalized Laguerre polynomials, 

{-Y v\x- v L^- v ix) = {-Yii\x-»L v -»{x) 
which, applied to the present case, allows us to write, when n q < m q 



(34) 



^•- W '[(M,/Rw fl ) 2 ] = (-)"«-»"« 



Mq 

fiu! n 



(rn q -n q ) 



xL^- n "[{M q /Tuj q f 



(35) 



From now on, we will always write the matrix element as displayed in_En.(3.14) above. Up to this point what we 
have is basically the generalization, for N vibration modes, of the resultsuCI for just one vibration coordinate. But the 
similarities between these and the results of Wingreen et al. in Ref. || lead us to believe that the TV-mode situation 
could somehow be reduced to a single vibration mode inelastic scattering event. This is in fact possible, if we impose 
the restriction of an Einstein band of phonons, that is to say oj q — luq for all modes q. In order to perform this 
reduction from N to just one "mode", we begin by considering only the last two phonon modes in Eq.(3.14) above, 
with the complete energy denominator: 



jv-rojv 
m N 



(9n) 



mjv-i— 



mjv— 



[(z - A N - 2 ) - {mN-xTiu)N-\ + m N Tiuj N )] 1 



(36) 



where we introduced the short notations z = e ki — e a + ^2 q (M^ /huj q ), A N _ 2 = mifkJi + m 2 Tiu 2 + • ■ ■ + m N _ 2 h UJ N-2 
and g q — (M q /huj q ) 2 . Having isolated these factors, we keep all summation indices m q fixed except for itin-i and 
toat, impose ujn-i = ujn = 1^0, and rewrite the double summation of these last two modes as 



(37) 



tun -1—0 mjv— 



which, by way of the general algebraic rule 



EE^ fe >' — E E ^p-m-p 

k=0 1=0 rn=Qp=a 



(38) 



can now be rewritten as 



G 



E 



J2Lr-^ p (9N-i)LZ_- m+p (g N ). 



(39) 



Appealing now to the sum rule (A7) for the generalized Laguerre polynomials, deduced in the Appendix below, the 
summation over the index p obeys the equality 



p=0 

TYl \9N-l+ 9N), 

and the contraction of the last two vibration modes summations originates the single summation 



E 

m=0 



i-r 



(z - A N - 2 ) - mTiujQ 



Lt N - 1+nN) - m (9N-i+9N). 



(40) 



(41) 



Keeping on contracting all Laguerre polynomials over the phonon modes, we end up with just one Laguerre polynomial 
and then the transition amplitude matrix clement can be written in a "one-mode" form: 



(l f ;ni,n 2 , ■ ■ ■ ,n N \T\ki;0i,0 2 ,- ■ ■ ,0 N ) = 
V lf , a V a , ki e-9 



n^) n ' /2 Ko- i/2 



E 



{ _ rL n m m {g) 

' [eki - e a + A - mhuj } ' 



(42) 



J2 q n q being the total number of final excited phonons into all the lattice vibration modes, g = ^2 a (M q /hoJo)' 2 



£ 9 3 9 andA = V K/^o). 



IV. SUMMING OVER FINAL PHONON STATES 



The transition probability will be given by the modulus square of amplitude above. Summing over all possible final 
phonon states (m, n 2 , • • • , njv| such that the total number of excited phonons is n, we have 



|(; /; n|T|fc i; 0)| 2 =\V lf , a \ 2 \V a , ki \ 2 e- 2 ° 



EE-E 



ni n 2 
oo 



n\\ n 2 \ 



E 



{ _ rL n- m(g) 



9_T 

n N \ 

2 



(43) 



However, in general x n /n\ = (—) n L~ n (x) and, taking advantage once more of the algebraic properties of the gen- 
eralized Laguerre polynomials, expressed in this instance by the sum rule (A10) in the Appendix below, we obtain, 
keeping in mind the restriction n\ + n 2 + ■ ■ ■ + raw = n, 

n ni n 2 n N 

^ ml n 2 l n N \ 

ni U2 tin 

= E E • • • E(-) niL «r (92)(-r 2 L-: 2 (92) 

™1 ™2 "JV 

■■■(-T n LZ n (9n) 



= (-) n L- n (g) = 



(44) 



Consequently, the transition probability from zero to all possible number n of final phonons (and from electron state 
ki to electron state If), imposing energy conservation in the overall process, will be given by 
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71 = 



]r|a /; n|T|fc i; 0)| 2 = 

oo n 

\Vi f , a \ 2 \V aM \ 2 e- 2s J2 ^(e fci - e,, - n^ ) 

{ _ )mL n- m{g) 



E 



efe, - e Q + A - m^wo] 



(45) 



for all the N phonon modes, but in a form that reproduces a one-mode situation, as sought. And if we now go from 
the domain of wavevectors ki and If, to the domain of incoming and outgoing electron energies, and e/, respectively, 
and taking matrix elements VJ /j(l and V a ,ki independent of the energy in the range of interest, we get a transmission 
matrix T(e,;, e/) given by 



E 

m=0 



[e< - e a + A - raftwo] 



(46) 



where = 27r^ fe . iV^feJ 2 ^ — efej and T/ = 2irJ2i f \V a ,i f \ 2 S(e — ei f ) are the partial widths of the intermediate 
resonant state due to the coupling to the continuum of electron states, initial and final. 
In this final result we may represent the generalized Laguerre polynomial as 



l (9) = E 



r(n + 1) 



(~9Y 



(47) 



and then, using the algebraic rule in Eq.(3.19), the factor inside the modulus square of Eq.(4.3) above can be rewritten 
as 



oo m 

EE 

m=0 j'=0 



(-)T(n + 1) 



(~9) J 



T(j — m + n + l)(m — j)\ j\ 



A — rnhwx 



o 



- ~ (_)(m 1+ m 2 ) r(Tl+1) ( )mi 



(48) 



However, if we take into consideration the presence of the poles of the gamma function T(— m 2 +71 + 1) at the points 
777-2 — ti + 1, 77 + 2, • • •, we conclude that the upper limit in the 7712 summation is, in fact, not 00 but 77. Consequently, 
we can write the transmission matrix above as (changing labels, from mi to m, and 7772 to j) 



T{e u £/) = TiT k e~ 29 2J —Sfe -e f - nhu Q ) 



71=0 



EE 

771=0 j=0 



(-)j'n! g m 
(n — ml 



e a + A - (m + j)huj Q } 1 



(49) 



which is the result expressed by Eq. (32) in Ref. || and by Eq. (10a) in Ref. [|. We may conclude, therefore, that 
the relationship between the end results of the two approaches goes beyond a numerical indistinguishability, and that 
a consistent description of inelastic resonant tunneling in heterostructures can also be achieved, without the need to 
use many-body Green's functions. 



V. CONCLUSIONS 



Our aim in the present paper was to clarify the possible validity of the antecipated connection between two different 
approaches to the same basic phenomenon, not only for one vibrational degree of freedom as done before, but for a 
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whole set of vibration modes, even if with the restriction of equal frequencies. Fundamental to the whole exposition 
were the algebraic properties of the generalized Laguerre polynomials. Furthermore, we believe that it may become 
possible, with this simpler method, to treat more general situations in a direct fashion, in particular the consideration 
of non-Lorentzian lineshapes, of non-linear electron-phonon coupling, and of more than one intermediate resonance. 
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APPENDIX: SUM RULES FOR THE GENERALIZED LAGUERRE POLYNOMIALS 



We give here proof of the sum rules between generalized Laguerre polynomials used above. A set of useful references 
regarding orthogonal polynomials is given below (Refs. ^l|-^7|). Starting from the generating function 

oo 

w(x, z) = (1 + z) a e- xz = c n {x)z n (Al) 

n=0 

expanded in powers of the complex variable z (such that \z\ < 1), the coefficients c n (x) are given by 

If w(x,z) 1 d n w{x,z) 

ler[(l + z) a e-*'] . f . 

— 1 2=0) ( AZ ) 



where the contour Y encloses the origin. Using Leibniz's rule for the n-th derivative of a product of two functions, we 
get 



Hence, we can write 



n\ ^ V r J (a-r)\ 

y> r(a + l) {-x) r 

^— ' Y(a — n + r + l)(n — r)\ r! 

r— v / \ / 

Lt- n) (x). (A3) 



w(x, z) = (1 + z) a e~ xz = L^- n \x)z n , (A4) 



using the definition of generalized (as opposed to "associated" ) Laguerre polynomials, where the upper index can be 
any complex number. 

Some convenient representations of the generalized Laguerre polynomials, which also remain valid for a a negative 
integer, are 



Y(a + r + 1) (n — r)\ r\ 



r=0 



1 J2 { -^(a + r + l) n _ r x r . (A5) 



n\ * — ' r 

r=0 



In this last expression we introduced the Pochhammer symbols, (a) r = Y(a + r)/Y(a). 

But now, with the help of the generating function above, we can prove the sum rule required. In fact, 
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{l + z) a e' xz ■ (l + zfe- yz 

OO OO 

= ££4 a - r) (*)^- s) 0/K + 

r=0 s=0 

•\a+/3„-(a:-H/)z 



= {l + z) a+l3 e 

OO 

= £ L( l a+/3 -") (a; + y)z n (A6) 



n=0 

Equating equal powers of the variable z, we conclude that r + s = n and consequently 



£ 4«~r ){x)L^r n \y) = L^-\x + y), (A7) 



r=0 



our first desired result. The summation index r can not attain values larger than n, since the lower index of a 
generalized Laguerre polynomial is always nonnegative. 

Next we prove the second sum rule used above. It is, in fact, a consequence of the first sum rule, for which the 
particular case a = (3 = gives 

n 

J24- r H Xl )L^ r n \x 2 ) = Li- n \x 1+ x 2 ). (A8) 

If we now consider only the first two factors in the product (4.2) above, we may write, for a fixed value of n (and of 
"3, • • ■ ,n N ), 

J2 L-^( gi )L-^(g 2 ) 

ni+n 2 =n-(n 3 -\ hnjv) 

n—(n3-\ h«jv) 

- 2^ L H WlJ L [ n -(„ 3 + . ,. +njv )_ ni ] (92) 

n 1= 

=^rn [ -^.;S )I (^+^)- (A») 

We repeat this contracting process over all N factors, ending up with just one final Laguerre polynomial, L~ n (g). In 
general, then, 

£ L-^( Xl )L-^(x 2 )...L-^(x N ) 

ni+n2H \-nN—n 

= L~ n (x), (A10) 
with x = xi + x 2 + • • • + xn and n = n\ + n 2 + ■ ■ ■ + tin- 
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